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Can Non-Gaussian Cosmological Models Explain the 
WMAP's High Optical Depth for Reionization? 
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ABSTRACT 

The first-year Wilkinson Microwave Anisotropy Probe data suggest a high optical 
depth for Thomson scattering of 0.17 ± 0.04, implying that the universe was reion- 
ized at an early epoch, z ~ 20. Such early reionization is likely to be caused by UV 
photons from first stars, but it appears that the observed high optical depth can be 
reconciled within the standard structure formation model only if star-formation in the 
early universe was extremely efficient. With normal star-formation efhciencies, cosmo- 
logical models with non-Gaussian density fluctuations may circumvent this conflict as 
high density peaks collapse at an earlier epoch than in models with Gaussian fluc- 
tuations. We study cosmic reionization in non-Gaussian models and explore to what 
extent, within available constraints, non-Gaussianities affect the reionization history. 
For mild non-Gaussian fluctuations at redshifts of 30 to 50, the increase in optical 
depth remains at a level of a few percent and appears unlikely to aid significantly 
in explaining the measured high optical depth. On the other hand, within available 
observational constraints, increasing the non-Gaussian nature of density fluctuations 
can easily reproduce the optical depth and may remain viable in underlying models 
of non-Gaussianity with a scale-dependence. 



1 INTRODUCTION 



The excess large scale polarization signal in the Wilkinson 
Microwave Anisotropy Probe (WMAP) data has allowed the 
first measurement of the optical depth to re-scattering by 
electrons with a value of r = 0.17±0.04 (Kogut et al. 2003). 
If interpreted as a sharp transition to a reionized universe 
from a neutral one, the redshift at which this transition hap- 
pens is = 17± 5 (Kogut et al 2003; Spergel et al 2003). 
Theoretical studies based on Cold Dark Matter (CDM) mod- 
els suggest that such a high reionization redshift — beyond 
what was previously suggested with Gunn-Peterson troughs 
in 2 ~ 6 quasars from the Sloan survey (e.g.. Fan et al. 2002) 
— generally requires significant star-formation early on, at 
z > 10 (Cen 2003; Haiman & Holder 2003; Sokasian et al. 
2003; Yoshida et al. 2003). 

While conventional models of reionization based on or- 
dinary stellar populations (Population II stars) formed in 
galaxies alone are still consistent with the measured opti- 
cal depth (Chiu, Fan & Ostriker 2003; Somerville & Livio 
2003), these models are not capable of explaining the high 
end of the allowed range for the optical depth. The best fit 
value of r = 0.17 is only achieved in models where reion- 
ization takes place with the maximal efficiency in terms of 



star-formation and associated ionizing photon production 
(Fukugita & Kawasaki 2003; Haiman & Holder 2003; Cen 
2003; Sokasian et al. 2003; Wyithe & Loeb 2003). Aherna- 
tively, with a combination of a high photon escape fraction 
and a top-heavy IMF in early galaxies, the resulting total op- 
tical depth can be increased to be marginally consistent with 
the WMAP data (e.g., Ciardi, Ferrara & White 2003). Al- 
though recent theoretical studies consistently indicate that 
the first stars formed out of a chemically pristine gas were 
rather massive (Abel, Bryan & Norman 2002; Bromm, Coppi 
& Larson 2002), hence supporting at least the latter condi- 
tion, it is unclear whether such conditions are realized in 
primeval galaxies over a sufficiently long period. It is there- 
fore intriguing to consider alternative mechanisms that do 
not require considerably high efficiency in star-formation or 
in ionizing photon production rate. Interestingly, such a sce- 
nario is possible even within the standard cosmological mod- 
els. For example, as we show later in the present paper, with 
erg equal to or larger than unity, the optical depth to reion- 
ization can reach within the range of the WMAP result since 
with a higher normalization, halos collapse slightly at earlier 
times and are more numerous at a given epoch when com- 
pared to models with a lower value for the normalization of 
the matter power spectrum. 
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In light of this, with currently favored values for cosmo- 
logical parameters, invoking non-Gaussianities in primordial 
density fluctuations may be an interesting possibility and de- 
serves further study. In the non-Gaussian models, one finds 
an increased fraction of high density peaks relative to Gaus- 
sian models, and these peaks are expected to collapse ear- 
lier than in the standard case (Peebles 1997; White 1998; 
Koyama, Soda & Taruya 1999; Robinson & Baker 2000; 
Willick 2000; Seto 2001; Mathis 2002). If star-formation 
is triggered in these collapsed halos, photons emitted from 
them may be able to, at least partially, reionize the universe 
very early on, causing a large total Thomson optical depth. 

Excess non-Gaussian fluctuations generally lead to for- 
mation of structures at late times and could overproduce the 
cluster abundance at low redshifts (Willick 2000). However, 
the constraint obtained from cluster abundance is rather 
uncertain, and previous studies based on models of pri- 
mordial non-Gaussianities (e.g., Verde et al. 2001) suggest 
that, in certain models where non-Gaussianities are asso- 
ciated with the density fleld, even a small primordial non- 
Gaussianity can signiflcantly affect the number counts of 
objects at high redshifts. Avelino & Liddle (2003) recently 
studied the evolution of the mass fraction in collapsed ha- 
los in non-Gaussian models and calculated a simple lower 
estimate for the redshift of reionization. 

In this Letter, we use a detailed model for the mass 
function in non- Gaussian models and compute the reioniza- 
tion history for some cases. To this end we compute the 
total Thomson optical depth as a function of redshift. We 
conclude that such non-Gaussian models, consistent with 
current constraints, can easily explain WMAP's results. We 
also argue that, while these models may eventually be con- 
strained with large scale structure observations, there is still 
some freedom to obtain the necessary non-Gaussianity with 
models which are strongly scale-dependent. 



2 CALCULATIONAL METHOD 

Following Haiman & Holder (2003) (hereafter HH), we study 
star formation and reionization using a simple toy model. 
Throughout the paper we assume the best fit flat ACDM 
model with power law fluctuations for the WMAP data: 
(n™,nA,n6,/i,cr8,n,) = (0.27,0.73,0.044,0.72,0.9,0.99). 
Given these cosmological parameters, we calculate the col- 
lapsed fraction of baryons using a spherical halo model. HH 
distinguished three different types of halos based essentially 
on the halo mass, which captures different nature of stars 
formed in them. The collapsed halo fractions for each halo 
type is 



-fcoll 



1 



dm mn{m, z) , 



(1) 



where mi and m^ take values based on the gas temperature 
within the halo such that 

m(.) = 0.94x10«Mq(A)"'(^)-''' 
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Figure 1. The fraction of mass in collapsed halos, between var- 
ious mass limits (see Section 2 for a summary), as a function of 
redshift for the standard Press-Schechter mass function and for 
the non-Gaussian mass function that includes a free parameter 
A that accounts for the departure from Gaussianity. Small values 
for A represent a significant non-Gaussianity while A > 100 is 
consistent with Gaussian fluctuations in the density field. 



where /i is the mean molecular weight: in the case of ionized 
gas with T > 10* K, fi — 0.6 and for neutral gas with 
r < 10* K, ^ = 1.2. The collapse fractions are divided to 
three halos between 100 K and 10* K (called Type II where 
molecular Hydrogen cooling dominates), between 10* K and 
2 X lO'^ K (Type la), and above 2 x 10^ K (Type lb). In both 
Type I halos, the dominant cooling mechanism to form stars 
involves atomic hydrogen lines. 

Using the collapsed mass fractions, we can now calcu- 
late the ionized HH filling factor, which is the same as the 
mean electron fraction Xe{z), by integrating for ionization 
produced by collapse at all higher redshifts, 

dz' Ui.^f^iz')+[l- Fuuiz')] X 



dz 



-{z') + eu 



dz 



Vmi{z',z). (3) 



The second term on the right hand side includes the fac- 
tor (1 — -Fhii) to allow for photo- ionization feedback, while 
Vmi[zon, z) is the volume per unit mass and unit efficiency 
ionized by redshift z by a single source that turned on at 
an earlier redshift of Zon > 2. In Eq. Q, = N-yf^fcsc 
represents the overall efficiency of star-formation and pho- 
ton emission, where /, is the fraction of baryons in the halo 
that turns into stars, A^'^. the mean number of ionizing pho- 
tons produced by an atom cycled through stars (averaged 
over initial mass function), /esc is the fraction of these pho- 
tons escape to IGM. In the following we use fiducial values 
of e for type la and lb halos as advocated by HH. We set 
en = 0, i.e. we do not consider star formation in "minihalos" 
in which the gas can cool via molecular hydrogen cooling. 
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It is straightforward to include Type II halos in our models, 
but the ionization history involving them is more compli- 
cated. For our purpose - test the effect of non-Gaussianity 
- it is unnecessary to introduce such additional complex- 
ity. Indeed, one aim is to see if early reionization could be 
obtained without invoking molecular hydrogen cooling. 

The evolution of the ionization front is calculated by 
solving the following differential equation 



at uh 



(4) 



where 0.3 ~ 2.6 x 10~^^ cm'' s~^ is the case B recombi- 
nation coefficient. Adopting a fixed value of as is equiv- 
alent to assuming a constant temperature of the ionized 
region. Strictly speaking, this is not true, and one needs 
consider the evolution of the thermal state of the ionized 
IGM. However, with atomic cooling the gas in the ionized, 
photo-heated regions quickly approaches this temperature, 
and as shown by HH this provides a good approximation. 
C = {njj) /fijj is the clumping factor, which we take to be 
at a constant value of 10. H{z) is the expansion rate, with 
//o[ilm(l + z)^ + ^Ia]^^^ , and N-y{t) is the rate of ionizing 
photon emission as a function of time. We take N^{t) to be 
a constant for t < 3.1 x 10^ years with a value of 3.7 x 10*® 
s~^ per solar mass of stellar content and decrease from this 
constant value as {t/S.l x lO^yr)"*'^ thereafter. We start 
integration at a sufficiently high redshift, and solve Eq. J^J 
up to Fhii = 1, at which point we consider reionization is 
completed and keep -Fhii fixed at 1 for redshift below this^. 
The optical depth is calculated by integrating the ionized 
volume fraction weighted by the number density such that 

t{z) = 1.08 J dz —caTnn{z)Fmi{z). (5) 

The factor 1.08 accounts approximately for the contribution 
of helium reionization (see HH) . 



3 NON-GAUSSIANITY 

The collapsed fraction can be calculated using a standard 
mass function, such as the Press-Schechter theory (PS; Press 
& Schechter 1974) based on Gaussian density fiuctuations, 
or the numerical simulated mass function by Jenkins et al. 
(2001), which was used by HH. To describe the linear power 
spectrum, we use the transfer function from Eisensetin & 
Hu (1998). We now modify the standard mass function by 
adding non-Gaussian fluctuations to the density field. To 
compare with various available constraints, we make use of 
a simple model for the non-Gaussian mass function which 
has been previously studied in the literature (Koyama, Soda 
& Taruya 1999; Willick 2000), though we note that, under 
specific models of primordial non-Gaussianity, more detailed 



^ Since this set of equations does not include evolution of IGM 
thermal state, it produce unphysical result (-Fhii > 1) if inte- 
grated past complete reionization, however the result is approxi- 
mately correct if the integration is limited to before the comple- 
tion of reionization. 



mass functions can be constructed (e.g., Matarrese, Verde & 
Jimenez 2000). We write the generalized form of the mass 
function as 



n{m, z) 



da{m, z) 



ma'^{m,z) dm 



a{m, z) 



(6) 



where the probability distribution function for density fluc- 
tuation, S, is written in terms of a dimensionless function 



Pi5\a) 



(7) 



where J ip{x)dx = 1, J xip{x)dx = and J x'^'ip{x)dx = 1. 
Here, cr(M, z) is the rms mass function on a scale corre- 
sponding to a mass m, and Sc is the critical density for col- 
lapse (see, e.g., Kitayama & Suto 1996 for details). The stan- 
dard Press-Schechter formula assumes a Gaussian probabil- 
ity distribution function (PDF) such that 'ip{x) = 2Ae~'' 
and A = {2'k)~^^^. To include a mild non-Gaussian tail to 
the Gaussian description, essentially a skewness, we consider 
a PDF which is based on the Poisson description, but modi- 
fied appropriately for the continuous variables involved (see, 
Willick 2000 for details). The function is written as 



Tp{x) 



: 2A 



r(\/Ax + A + 1) 



(8) 



such that A, which is a function of redshift, is a free- variable 
that captures the non-Gaussian nature of the PDF; As A is 
increased, the mass function approaches to the PS value for 
Gaussian distribution of density fluctuations. A is chosen to 
return the correct normalization: J mn{m, z)dm — p{z). 

Note that in our non-Gaussian mass function, skewness, 
5*3, scales as A~^, which explains why non-Gaussianity in- 
creases with decreeising A. In terms of the moments of the 
primordial density field, one can write skewness as 



Ssim, z) 



jJL3{m, z) 
a'^ijn, z) 



(9) 



where ^3 is the third moment of the density perturbations, 
scaled to a redshift of z, when e is the parameter related 
to the quadratic correction (Matarrese, Verde & Jimenez 
2000). If we ignore for now the mass-, or equivalently scale- 
dependence, of primordial non-Gaussianity, then /is ~ 6(t'*, 
and assuming linear growth of fiuctuations, Sz{z) ~ (ye/G(z) 
where G{z) is the growth factor. The latter scaling can be 
understood based on the fact that p.i(m,z) and a^{m,z) 
scale with redshift as G^(z) and G^{z), respectively. Since 
in our mass function, 5*3 ~ A^^, we expect A to behave as 
~ /6)G{z) as a function of redshift. 

In Fig. 1, we plot the collapsed fraction of these halos for 
the standard PS mass function and the non-Gaussian mass 
functions with A = 10 and A = 1 as a function of redshift. 
We note that these values could be a function of redshift, 
although we have taken A to be a constant, for purposes of 
the present discussion. As A is decreased, or models become 
strongly non-Gaussian, the fraction of mass in collapsed ha- 
los increases at a fixed redshift. This illustrates the basic idea 
behind our calculation that with increasing non-Gaussianity 
of the density field, there are more high overdensity peaks 
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Figure 2. Left: The volume ionized fraction, Xe{z) = Fiiii{z) as a function of redsliift. Middle: The Thomson optical depth for electron 
scattering. The models are based on "ordinary" star-formation in Type I halos only. The Press-Schechter mass function leads to an optical 
depth with a value of 0.097 while the non-Gaussian mass functions lead to values of 0.11 and 0.13 for A = 10 and A = 1, respectively. 
The two black dashed lines in the middle and right panels indicate the range of the first year WMAP result. Right: Same as the middle 
plot, but with (Tg = 1 and keeping all other parameters the same. 



which collapse earlier. These early collapsing halos can form 
stars at an earlier time than in the standard description and 
reionize the universe such that the total optical depth is 
higher than previously estimated in the Gaussian models. 

In Fig. 2, we show the ionization fraction and optical 
depth for the models with the PS mass function and with 
the non-Gaussian mass function with A = 10 and A = 1. In 
addition to our standard cosmological values, to illustrate 
the increase in the optical depth by halos collapsing at an 
earlier time, we also calculate the optical depth in a model 
where ag = 1 (Fig. 2, right plot). We emphasize that such 
high values of ug easily make the optical depth large even 
within the Gaussian model. 

For a given set of cosmological parameters, higher op- 
tical depths can be achieved with the non-Gaussian mass 
functions. As shown in Fig. 2, compared with the standard 
Gaussian PS function, there is only a moderate increase 
in optical depth for the A = 10 model (about 0.01 in r). 
This difference is in fact comparable in size with the differ- 
ence between the PS and Jenkins mass functions: it is well 
known that compared with the PS mass function, the Jenk- 
ins et al. (2001) function predicts an excess of halos at the 
upper end of the halo mass. Our calculated optical depth 
with the Jenkins et al. (2001) mass function of ~ 0.13 is 
consistent with previous calculations on the optical depth 
based on Pop II stars alone (e.g., HH; Fukugita & Kawasaki 
2003). To increase the optical depth beyond this value, one 
can either add more Population III stars at earlier epochs 
(Cen 2003; Wythie & Loeb 2003; Sokasian et al. 2003) or, as 
we have considered, by considering smaller A value (greater 
non-Gaussianity) while resorting only ordinary stellar pop- 
ulations. 

A significantly non-Gaussian model with A ~ 1 in- 
creases the optical depth by a large fraction (about 0.04 in 
r). From Eq. |5]and the expected scale with redshift, such a 



non-Gaussianity implies a primordial non-Gaussianity with 
e ~ 5 X 10~^. Non-Gaussianities at this level can soon 
be tested with observations of galaxy counts at redshifts 
of 5, though low redshift large scale structure observations 
will not be helpful for this purpose (Verde et al. 2001). At 
present, according to Willick (2000) (see his Fig. 11), for 
f2m = 0.27, the low redshift galaxy cluster abundance is 
consistent at the 90% confidence level with a non-Gaussian 
mass function of A < 14. The cluster abundance at low- 
redshift is not a strong discriminator of non-Gaussian mod- 
els, as argued by Verde et al. (2001). Note that only a con- 
straint on the upper limit of A is available from the cluster 
abundance in Willick (2000), though in Koyama et al. (1999; 
also, Robinson, Gawiser & Silk 1998) A ~ 3 is preferred for 
ACDM cosmological models. Again this is not a strong con- 
straint of non-Gaussianity at redshifts of order 30. 

Apart from low and moderate redshifts observations, 
non-Gaussianity can also be constrained, at a redshift of 
1100, from observations of cosmic microwave background 
anisotropics (e.g., Komatsu et al. 2003). Current constraints 
indicate a limit of 134 on /ni, at the 95% confidence level; 
this parameter is related to the quadratic corrections to po- 
tential field and is different from e we have been using to 
describe primordial non-Gaussianity based on corrections to 
the density field directly. Through the Poisson function, one 
can relate the two parameters such that at large scales of or- 
der 100 Mpc, roughly, /ni ^ 2 x lO'^e, and the upper limit on 
/ni of 134 from WMAP corresponds to e < 7 x 10~^ and al- 
lows for significant non-Gaussianity, rendering a lower limit 
on A at a redshift ~ 30 of 0.08. While cluster counts from 
low redshifts have constrained the upper limit and CMB 
observations have constrained the lower limit, the allowed 
range is still significant enough to produce substantial non- 
Gaussian fluctuations and to easily reproduce the observed 
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optical depth with early formation of structures and first 
stars. 

Note that with completed WMAP results and with 
Planck, the limits on non-Gaussianity will be improved and, 
naively, small values for A at 2 ~ 30 may be ruled out. Such 
a limit, however, cannot fully constrain non-Gaussianities, 
since one might still circumvent this by considering scale- 
dependent non-Gaussian models. For example, one can ap- 
ply the significant non-Gaussianity to scales which are re- 
lated to star-formation at redshifts above 10 and still allow 
such non-Gaussianity to remain today. It will then affect the 
number counts of ~ 10* Mq halos at low redshifts, but CMB 
anisotropics at large scales will remain nearly Gaussian. The 
same models can also be tuned such that cluster counts at 
low redshifts remain Gaussian also. We note that it is not 
difficult to envision how such models can be constructed 
with examples including non-standard inflationary models 
(Wang & Kamionkowski 2000; Bartolo, Matarrese & Riotto 
2002; Bernardeau & Uzan2003) and models with primor- 
dial black hole formation proposed by Bullock & Primack 
(1997). While constructing a concrete model for producing 
such scale-dependent non- Gaussian primordial fluctuations 
is beyond the scope of the present study, we argue that these 
models can also be constrained eventually. The best way 
to probe the level of non-Gaussianity suggested above, and 
at scales related to this problem, may be to study number 
counts of galaxies at the highest redshifts possible (Verde et 
al. 2001). Another possibility involves a measurement of the 
non-Gaussianity in the surface brightness fluctuations pro- 
duced by flrst stars and proto-galaxies in the near-infrared 
background and potentially detectable with upcoming deep 
observations from Near-IR Imaging Camera on the ASTRO- 
F mission (Cooray et al. 2003). Studying these constraints on 
the presence of non-Gaussianities at highest redshifts related 
to formation of flrst objects is clearly needed to understand 
the extent to which the reionization process is affected. 



4 SUMMARY 

We have examined how the observed high optical depth 
is achieved in cosmological models with non-Gaussian pri- 
mordial density fluctuations. In these models, high over- 
density peaks due to the non-Gaussianity collapse early on 
and may harbor first stars at redshifts above what is gen- 
erally expected in the standard Gaussian model. For mild 
non- Gaussian fluctuations at redshifts of 30, the increase in 
optical depth is at the level of 0.01, and it is unlikely to ex- 
plain the observed value of optical depth by this mechanism 
alone. To obtain significant increase in the optical depth, 
the first collapsing halos associated with the primordial non- 
Gaussianity must contribute to reionization significantly. A 
substantial volume fraction of the IGM may be reionized if a 
number of Population III stars are formed in these early ha- 
los. We found that increasing the degree of non-Gaussianities 
can easily reproduce the high optical depth. Such models 
are still viable within constraints from both cluster number 
counts at low redshifts and CMB anisotropics at a redshift 
of 1100. While constraints from these two redshift ends will 
improve, models with scale-dependent non-Gaussianity are 



possible to explain the optical depth. Such models will only 
be constrained by probes of non-Gaussianity related to mass 
scales of 10* Msun at redshifts of 30 and will involve, for ex- 
ample, study of first galaxy counts at very high redshifts. 

The ongoing operation of WMAP will pin down a 
presice value for the total optical depth. For a longer-term, 
post- WMAP CMB polarazation experiments such as Planck 
will probe the reionization history (Kaplinghat et al. 2003), 
and detection of the second-order polarization anisotropics 
can place additional constraints on details of the reioniza- 
tion process (Liu et al. 2001). Data from these future ob- 
servations will improve our picture of reionization and will 
enable us to further distinguish theoretical models related 
to the first generation of luminous structure formation. 



ACKNOWLEDGMENTS 

We thank organizers of the US-Japan Workshop on the SZ 
effect, where this paper was conceived, for a wonderful con- 
ference. This work is supported by the Sherman Fairchild 
foundation and DOE DE-FG 03-92-ER40701 (at Cahech) 
and by NSF PHY99-07949 (at KITP). NY acknowledges 
support from Japanese Society of Promotion of Science Spe- 
cial Research Fellowship (Grant 2674). NS is supported by 
Japanese Grant-in-Aid for Science Research Fund of the 
Ministry of Education, No. 14340290. 



REFERENCES 

Abel, T., Bryan, G.L., & Norman, M.L., 2002, Science, 295, 93 
Bartolo, N., Matarrese, S., Riotto, A., 2002, Physical Review D 
65, 103505 

Bernardeau, F. & Uzan, J. P., 2003, Physical Review D 67, 121301 
Bromm, V., Coppi, P.S., & Larson, R.B., 2002, ApJ, 564, 23 
Bullock, J. S. k. Primack, J. R., 1997, PRD, 55, 7423 
Cen, R., 2003, ApJ, 591, 12 

Ciardi, B., Ferrara, A., White, S.D.M., 2003, submitted to MN- 
RAS 

Cooray, A., Bock, J. J., Keating, B., Lange, A. & Matsumoto, T. 

2003, ApJ submitted, astro-ph/0308407 
Eisenstein, D. & Hu, W, 1998, ApJ, 496, 605 
Fan, X., Narayanan, V. K., Strauss, M. A., White, R. L., Becker, 

R. H.., Pentericci, L. & Rix, H. 2002, AJ, 123, 1247 
Fukugita, M. & Kawasaki, M. 2003, MNRAS, 343, L25 
Haiman, Z., Abel, T., & Madau, P. 2001, ApJ, 551, 559 
Jenkins, A., Frenck, C. S., White, S. D. M., et al. 2001, MNRAS, 

321, 372 

Kaplinghat, M., Chu, M., Haiman, Z., Holder, G. P., Knox, L., 

2003, ApJ, 583, 24 
Kitayama, T. & Suto, Y. 1996, ApJ, 480, 493 
Komatsu, E., et al, 2003, ApJS, 148, 119 
Kogut, A., et al, 2003, ApJS, 148, 161 

Liu, G.-C, Sugiyama, N., Benson, A. J., Lacey, C. G., Nusser, 

A., 2001, ApJ, 504, 516 
Mathis, H. 2002, PhD thesis, University of Toulouse 
Matarrese, S., Verde, L. & Jimenez, R. 2000, ApJ, 541, 10 
Koyama, K., Soda, J. & Taruya, A. 1999, MNRAS, 310, 1111 
Peebles, P. J. E. 1997, ApJ, 483, LI 
Press, W. H., Schechter, P. 1974, ApJ, 187, 425 [PS] 
Robinson, J. & Baker, J. 2000, MNRAS. 311, 781 
Robinson, J., Gawiser, E., & Silk, J. 2000, ApJ. 532, 1 



© 0000 RAS, MNRAS 000, 000-000 



6 Chen et al 



Seto, N., 2001, ApJ, 553, 488 

Shapiro, P. R. & Giroux, M. N., 1987, ApJL, 321, 107 
Sokasian, A., Yoshida, N., Abel, T., Hernquist, L. Springcl, V. 

2003, astro-ph/0307451 

Somerville, R. S. & Livio, M. 2003, |astro-ph/0303017| 

Spergel, D.N., et al, 2003, ApJS, 175 

Wang, L. & Kamionkowski, M., 2000, PRD 61, 063504 

Willick, J. 2000, ApJ, 530, 80 

Wyithe, S., & Loeb, A., 2003, ApJL, 586, 693 

Verde, L., Jimenez, R., Kamionkowski, M., & Matarrese, S. 2001, 

MNRAS, 325, 412 
Yoshida, N., Abel, T., Hernquist, L., & Sugiyama, N., 2003, ApJ, 

592, 645 



© 0000 RAS, MNRAS 000, 000-000 



